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Abstract 

Necessary conditions for a field theoretic equation of motion to 
be the consequence of variation of an infinite number of inequivalent 
Lagrangians are examined. 







1 Introduction 



In this article, an avenue of investigation suggested by the Universal Field 
Equations proposed in U is explored. Two of the striking properties of these 
equations are that they possess an infinite number of inequivalent Lagrangian 
formulations, and are completely integrable @, a feature which was antici- 
pated by the presence of an infinite number of conservation laws implied by 
the Lagrangian property. If, as is the case here the Lagrangian does not de- 
pend explicitly upon the field, the Euler variational equation takes the form 
of a conservation law. Also there exist many different conserved stress ten- 
sors. These attributes suggest that they can be considered as a halfway-house 
between ordinary Lagrangian and Topological Field Theories. The first two 
sections are devoted to an analysis of certain circumstances under which an 
equation of motion for a scalar field, containing no derivatives higher than the 
second in space-time of dimension two or greater admits an infinite number 
of inequivalent Lagrangians. (Two Lagrangians which differ by a divergence 
or are constant multiples of each other are of course equivalent.) Hundreds 
of papers have been devoted to this topic, the inverse problem in the calculus 
of variations, in the case where only time dependence is considered [[J , but 
comparatively little has been done in the case of field theories. It will be 
demonstrated that while not generic, the situation where the inverse prob- 
lem has a non-unique solution is not uncommon. However, if the additional 
requirement of Lorentz or Euclidean invariance is imposed, the possibilities 
are rather limited, and this rules out most such equations for any possible ap- 
plication to physics. The final section establishes integrability by the method 
of Legendre Transforms. 

The Universal Field Equation 

detfj^ =0, i = l-d. (1.1) 

\ dx j dxjdxk ' 

in d dimensional spacetime of which the Bateman equation 

( d<j) 2 d 2 (p { d(j) 2 d 2 (j) d(j)d(j) d 2 cj) 

[ dx' dt 2 [ dt' dx 2 dxdt dxdt ' 1 ' ] 

is the two-dimensional case admits an infinite number of Lagrangians, the 
construction of which is given in 0. (for the Bateman case, any Lagrangian 
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density of the form 



where F is an arbitrary function will do)! 

We shall argue that a necessary condition that the following generalisation 

of o); 

1 0. (1.4) 



E^(- 



-dx k dxi J v 

where G is a d x c/ matrix with elements functions of first derivatives J^- and 

g^ 2 ^. the Hessian matrix, regarded as an equation of motion, possesses an 
infinite number of inequivalent Lagrangian formulations is that det G is of 
rank at most d — 1. In the case of (|1.1|) the relevant matrix is of rank one, and 
this is overkill! We shall demonstrate this, first by examining the situation 
for d = 2 and d = 3, then by invoking a linearisation by means of a Legendre 
transform suggested by 0, infer the general result. 



2 Infinitely many Lagrangians in 2 dimen- 
sions. 

Consider a scalar field theory with Lagrangian £(&), i = 1,2, where <pi 
denotes J^-. The equation of motion is 

i=2 3=2 Q2p 

E E nj-r^ = Gu&i " 2GW12 + G 22( /> 22 = 0, (2.1) 

i=lj=l Cl( Pi ( Pj 

in the notation of (|1.4| ). From this formulation it is easy to make the identi- 
fications 

d 2 C d 2 L d 2 C 

7^2 = A (^)Gn; j-— = -\(<f>i)G 12 ; = A(^)G 22 . (2.2) 

o<K #102 dcpi 

where a proportionality function A has been introduced which will be deter- 
mined up to a constant factor if £ is unique. Now if the equation of motion 
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Q2.1D is given, there are two integrability constraints arising from ( f2.2| ) of the 
form 

902 »01 O02 901 

§G 12 + ^G 22 + A(^2 + ^2) = 0, (2.4) 

902 901 902 901 

to be satisfied. In general, if these equations are linearly independent they 
can be integrated to determine A uniquely in terms of the functions and 
their derivatives. If, instead, they are linearly dependent, then A retains a 
certain arbitrariness, which results in an arbitrariness in the construction 
of C from (|2.2|) . Necessary and sufficient conditions for this eventuality are 
easily obtained as 
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det[^ L n )= Q ( 2 - 5 ) 
0. (2.6) 



90i V Gn 

These are evidently satisfied by the Bateman equation ( |1.2| ) for which 

Gn = 4>l Gu = 0i0 2 ; G 22 = 4>l (2.7) 

More general solutions appear to be difficult to obtain, and in what follows 
we shall restrict ourselves to the case where gij is quadratic in 0^. For the 
choice ( |2.7D the equation to determine A ( |2.4| )is simply 

0i^T + 2 ^ + 3A = O. (2.8) 

901 902 

This simply means that A is a homogeneous function of 0j of degree -3, a 
result which may be verified to be in accord with ( |1.3| ). 



3 Infinitely many Lagrangians in 3 and more 
dimensions. 

The starting assumption of the previous section must be modified to obtain 
a nontrivial extension of these results to 3 dimensions. If the equation of 
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motion is simply linear in the second derivatives, with coefficients functions 
of the first derivatives of the fields, then the constraints analogous to (|2.4j ) 
are overdetermined. Instead we start with a Lagrangian density of the form 

£ = EEi ij W^i (3.i) 

i=l j=l 

Here the coefficients L y ' depend only upon the first derivatives of <p so 
that £ depends only linearly upon the second derivatives of the field. This 
is a redundant description, as three of these coefficients may be removed by 
removing divergences to eliminate equivalent Lagrangians. The equation of 
motion takes the form 

d 2 L ij 

anr^^H ~ Mik) = °- ( 3 - 2 ) 
i,j,k,i d( Pk<Pi 

Now compare this equation with the equation of motion, assumed given in 
the form 

^ y] £ipq£jrsG((j)k)ij4>pr<f)qs (3-3) 
i,P,qj,r,s 

Compatibility of those two equations requires the existence of a proportion- 
ality factor A such that 

d 2 L pr 

Y Y e ipQ € jrs o t 7 = ^Gi.j (3-4) 

p,q r,s Vtyqtys 

In this formalism it is easy to see that there are three equations of consistency 
guaranteeing the integrability of ( |3.4| ), namely 

gA(A Gy )=0. (3.5) 

The argument proceeds along the same lines as before; If the equations ( |3.5| ) 
admit a solution for A which is not unique (discounting a trivial proportion- 
ality constant) then a necessary condition is that det G = 0, i.e. the matrix 
of coefficients Gij must have rank 2 at most. There is also a differential equa- 
tion to be satisfied. At this point we introduce the simplifying assumption 
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that the G's are quadratic in derivatives in order to exhibit two solutions. 
The first 

Gij = (j)^ (3.6) 

recovers the Universal Field Equation of [Q for the case of 3 dimensions. The 
matrix G is here of rank 1 , a stronger requirement than necessary. The gen- 
eral solution of (|3.5| ) determines A to be an arbitrary homogeneous function 
of 4>i of weight - 4. 

The second is constructed using the following matrix representation for 

G; 

( 4>l + 4>l faifa + 2 ) 2 (0i + 3 )\ 
G= U<t>i + 4>2) + M<h + fo) • (3.7) 
\M<Pi + h) Mfa + fo) 4>l + (t>l 1 

In this case G is of rank 2, and the coefficients are such that the equations 
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( |3.5|) are equivalent to two independent equations; 

+ A = 0, (3.8) 

+ A = 0. (3.9) 

These equations admit the symmetric solution, among many possible solu- 
tions, 

A = {(ft + <Pl + <t>\ - 0102 - 0203 - 030l) 1+fe (01 + 02 + fof \ (3.10) 

which depends upon a single parameter k. The Lagrangian £ is then de- 
termined up to divergences in terms of the three functions L 11 , L 22 , L 33 
(setting the others to zero). For each choice of k the equations to for these 
functions may be solved, giving rise to an inequivalent set of Lagrangian 
densities parametrised by k. Two specific cases are diplayed in an appendix 
to this paper. 

This example may be used to deduce an automatic construction of equa- 
tions of motion in d dimensions which admit an infinite number of distinct 
Lagrangian formulations. 

Consider the following equation of motion; 

E^ Ad j(/^).. = o, (3.H) 

K 1 
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where 

Adjf ° 2 ^ ) -dctf ^ )( 9 ' ^ Y l (3 12) 

^ dxkdxi ' v ^ dxkdxi ' ^ dx^dxi ' v 

is the adjugate matrix, and is constructed in the following manner. 

Take (f)[ = P((4>i) to be a permutation of the components of the vector 

02; • • • 4>d} by an element P of the permutation group in d dimensions. 

Then take G^- to be a linear sum with arbitrary constant coefficients of terms 

of the form <f)f cpj, taken over subset of not more than d — 1 elements P of 

the permutation group. This will guarantee that det G is of rank at most 

d — 1. The case which corresponds to the three dimensional situation above 

is obtained by taking a sum over cyclic permutations. Take a differential 

operator representation for Pj as 

d<p j+1 d(pj + 2 dtpj 
where the indices are defined in an obvious cyclic manner and take 

Gij = Pi(0i)Pi(0i) + P 2 (^)P 2 ((0i) + • • • + Pd-ii^Pd-Mj). (3.14) 

Then, introducing A in an obvious generalisation of the 3 dimensional situa- 
tion above the d—1 consistency equations to determine A become 

P n \ + A = 0, n = 1, 2, . . . , d - I. (3.15) 

Suppose ujj is a <fth complex root of unity and we define u> = 1. Then the 
function fj = Z)f =1 u;]~ fa is an eigenfunction of the differential operator Pi 
with eigenvalue uj d ~ x = u*. Then a solution to the equations ( |3.15| ) dependent 
upon an arbitrary parameter k is evidently given by 

A = J II VVo (3-16) 
in analogy with (|3.10|). If d is composite, then there is an obvious way to 



include additional parameters in this type of solution, by multiplication by 
arbitrary powers of other zero eigenfunctions. For each such A, the equations 
to determine a Lagrangian will be integrable, and thus an infinite number 
of Lagrangians, parametrised by k will exist. Equations of this sort, unlike 
the Universal Equation, have lost Lorentz invariance and retain only discrete 
transformations among the space-time variables. This means that as far as 
physics is concerned, such equations can be discounted. 



6 



4 Legendre Transforms 



Further insight into the circumstances under which an infinite number of 
Lagrangian descriptions of a theory exist is afforded by the Legendre Trans- 
form, which was used in [[| to linearize the Universal Field Equation. This 
transform, which is clearly involutive has the flavour of a twistor transform, 
and also works for equations of the type ( |3.11| ). The multivariable version of 
this transform runs as follows [|J Introduce a dual space with co-ordinates 
£i, % — 1, . . . , d and a function Q(£i) defined as follows 

<f>(xi,x 2 , ...,x d ) + Qfa, f 2 , ...,&) =xi£x + X2C2 + x d i d . (4.1) 

C ^ 9 Q W (A 0\ 

6= a£' X * = W v " (4 ' 2) 

To evaluate the second derivatives faj in terms of derivatives of Q it is conve- 
nient to introduce two Hessian matrices; $, Q with matrix elements ^ and 
= gjj respectively. Then assuming that $ is invertible, $Q = 1 and 

^ "IS-*)., ^ = (•"%. (4-3) 



dxidxj <K,i)^j 

The effect of the Legendre transformation upon the equation (|3.11| ) is 
immediate; in the new variables the equation becomes simply 

E^HpO, (4.4, 

a linear second order equation for q. In the case discussed above, with a 
permutation operator representation for Pj now as 

p ' = ^ + ^ + - +( *i (4 ' 5) 



then equation (4.4) becomes 



E^E (PiP s -Piu)Q = 0. (4.6) 
i=i j=i 
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Now the crucial point is that this equation does not affect the dependence 
of the function Q upon the particular combination represented by a zero 
eigenfunction of all Pj, i.e. 

j=d—l i=d 

n (E<4 -1 &) ( 4 - 7 ) 

j=0 i=l 

. This means that it is effectively an equation in a manifold of dimension 
d — 1 instead of d. Thus the phenomenon of the existence of an infinite 
number of Lagrangians in the original space, is related to the fact that the 
linear differential operator in the transform space acts only upon a subspace 
of the transform variables . 



5 Appendix: Examples 

Two examples of Lagrangians for the equation ( |3.3| ) £1X6 clS follows: Corre- 
sponding to parameter k = — 1 

3 

C = - Mi(<l>i® l °9($) - 0i+i0i+ 2 ) (5.1) 
i=i 

where 

$ = 01 + 02 + 03 

and the indices are defined cyclically. The Lagrangian corresponding to pa- 
rameter k = is given by 



£ = + 03 - 20?) + 90?(0 2 + 3 ) + 901(0* + 

+ 450 2 0!0 3 (0^ + cf>l) - 15(0*0? + 4^ + 0^04 + 0^4) 

+ 30020^03(3(0201 - 0203 + 0103) ~ Hi + 0s))) (5-2) 

+ cyclic replacements. 
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